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Abstract 

The optical potential of halo and weakly bound nuclei has a long 
range part due to the coupling to breakup that damps the elastic 
scattering angular distributions. In order to describe correctly the 
breakup channel in the case of scattering on a heavy target, core recoil 
effects have to be taken into account. We show here that core recoil 
and nuclear breakup of the valence nucleon can be consistently taken 
into account. A microscopic absorptive potential is obtained within 
a semiclassical approach and its characteristics can be understood in 
terms of the properties of the halo wave function and of the reaction 
mechanism. Results for the case of medium to high energy reactions 
are presented. 
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1 Introduction 



Since the advent of Rare Isotope Beams (RIBs) P elastic nucleus-nucleus 
scattering with a radioactive projectile |2I is a reaction which has been studied 
to a large extent in the attempt to find characteristics that would be typical 
for a weakly bound nucleus and would help understanding new phenomena 
such as the existence of the halos 0. It has been established that for halo 
projectiles breakup of the valence particles is responsible for a damping in 
the elastic angular distribution starting from around 5° at medium to high 
energies. 

All theoretical methods used to describe the above mentioned reactions, 
require at some stage of the calculation the knowledge of the nucleus-nucleus 
optical potential. The optical potential is the basic ingredient for the de- 
scription of elastic scattering, but it is important also in transfer and breakup 
calculations, since one needs to take into account the core quasi-elastic scat- 
tering by the target while the valence neutrons are transferred or breakup. 
For example in the case of some two-neutron halo nuclei such as ^^Be or ^^Li, 
their cores of ^^Be and ^Li are themselves weakly bound nuclei. For multi- 
nucleon transfer reactions planned in order to obtain heavy exotic nuclei it 
will be important to have the appropriate optical potentials which include 
breakup and which should be used in the intermediate steps of the reaction. 
In the charge exchange reaction ^^B(^LiJ BeY^Be the halo nucleus-nucleus 
optical potential necessary to describe the final channel ^ has a volume 
part obtained with a double folding plus a very diffuse surface term fitted 
phenomenologically to reproduce the final channel angular distribution. The 
effect of the surface potential reduced the absolute cross sections by about 
50%jlj, in accordance with the experimental data and it can be interpreted 
as due to the halo breakup. 

Microscopic optical potentials for a halo projectile have already been stud- 
ied by many authors, and a review of the present situation can be found in 
13 in] . One of these methods consists in starting from a phenomenologically 
determined core-target potential and then to add the effect of the breakup 
of the halo neutron. This process leads to adding a surface part to the core- 
target potential. This new surface peaked optical potential has been seen 
to have a quite long range which reflects the properties of the long tail of 
the halo neutron wave function. Such kind of potentials are often called 
dynamical polarization potentials [Uj-fTTj. 

In a recent contribution we proposed a new approach [H] to the calculation 
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of the imaginary part of the optical potential due to nuclear breakup. It is 
based on a semiclassical method described by Broglia and Winther in [THI IT^ 
and used also by Brink and collaborators fK]\ 1^ to calculate the surface 
optical potential due to transfer and on the Bonaccorso and Brink model for 
transfer to the continuum reactions [221-1211 j the idea being that breakup 
is a reaction following the same dynamics as transfer but leading mainly to 
continuum final states when the incident energies per nucleon are higher than 
the average nucleon binding energy. The calculations were almost completely 
analytical and a simple, approximated formula was obtained which helped 
understanding the origin of the long range nature of the potential and its 
dependence on the incident energy as well as on the initial neutron binding 
energy. The characteristics of our potential were consistent with those of 
potentials obtained with other methods, in particular the eikonal method of 
Canto et al.[II)| and application to the description of experimental data were 
encouraging .5j. 

However it is very well known that for heavy targets recoil effects which 
give rise to the so called Coulomb breakup, are important and actually dom- 
inant for a neutron halo [23 1211 ISOl IHI] in the breakup cross section and 
therefore one wonders what would be their effect on the elastic scattering. 
Some work has already been published in order to calculate the optical poten- 
tial due to Coulomb breakup at low [EI, [IS] or intermediate energies [2Z1I2H]- 
However paper [2H1 deals with proton halo breakup which, as it has been de- 
mostrated in Ref. [21] has to be treated with great care when compared to 
neutron breakup. Therefore the methods used in |28j to include breakup in 
the optical potential are not expected to be applicable to the neutron breakup 
case. We will show in this paper that the method used in Ref . [S] can be used 
in the case of Coulomb breakup as well and that the corresponding formal- 
ism, which is appropriate to reactions performed at medium to high energies, 
well above the Coulomb barrier, is consistent with and joins continuously to 
the formalism used by other authors ^3] at lower energies. 

2 Theory 

The method we use here is the same as in [3] and it is based on the extraction 
of an optical potential from the calculation of a phase shift. 

The elastic scattering probability is Pei = IS'atatP, given in terms of the 
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nucleus-nucleus S-matrix. We know that 

\SMb)\' = e-'^'^'l (1) 

In a semiclassical approximation 18j, the imaginary part of the nucleus- 
nucleus phase shift 6i is related to the imaginary part of the optical potential 
by 

m = -- {Wy{v{t)) + wsm)) dt (2) 

where the volume potential is responsible for the usual inelastic core-target 
interaction, while the surface term takes care of the peripheral reactions like 
transfer and breakup. r(t) = he + vt is the classical trajectory of relative 
motion for the nucleus-nucleus collision. 

According to [5J the surface optical potential Ws(r(t)) due to breakup 
can be related to the breakup probability by 

/ + 00 fi 
Ws{r{t))dt = --P,M (3) 
-oo Z 

where P^^^ = J2iPi are the breakup probabilities in the various channels i. In 
order to obtain the surface imaginary potential Eq. (jS)) should be calculated 
as an identity in the distance of closest approach, which amounts to require 
that Ws{r) be a local, angular momentum independent function. We remind 
the reader that since we are using a semiclassical method, the non locality, 
which is in principle a characteristic of microscopic optical potentials has 
been transformed into an energy dependence |35j . 

Eq.Q can also be derived in a straightforward way from the time de- 
pendent scattering Schrodinger equation for the elastic channel probability 
density function in presence of a complex potential E2] • In the tradi- 
tional formulation the index (i) stands for stripping and pickup to bound 
states and in Ref. we extended it to hold for breakup reactions in which 
the final neutron state is in the continuum. Nuclear breakup of both absorp- 
tive and diffractive type was included and here we will include also Coulomb 
breakup. The justification of the use of Eq.© to calculate the imaginary 
potential due to nuclear breakup was simply given by the analogy between 
breakup and transfer as expressed by the transfer to the continuum model 
introduced in Refs. [221-1211 • There it was shown that the formalism for trans- 
fer to bound states goes over transfer to the continuum in a natural way if 
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the kinematics of the reaction is taken into account correctly within a time 
dependent approach which ensures neutron energy conservation. 

On the other hand Brogha and Winter in ^Hl d| pointed out the fact 
that the same formahsm could be extended to include core recoil. In the case 
of breakup of weakly bound nuclei we have shown in Ref . Oil that core 
recoil is responsible for the Coulomb breakup and that nuclear and Coulomb 
breakup give rise to negligible interference effects. Therefore we argue here 
that the imaginary part of the optical potential due to Coulomb breakup can 
also be calculated by Eq.Q where now one of the pi probabilities will be 
that of Coulomb breakup of the valence nucleon. 

Then, using Eq.Q and ©, in (P) the nucleus- nucleus S-matrix, in the 
case of a halo projectile, can be written as 

\Snn\' = l^cTre-^^- (4) 

where Sct takes into account all core-target interactions while the term 
Q-Pbup depends only on the halo neutron breakup probability. For a halo 
nucleus at high incident energy the transfer probability is going to be much 
smaller than the breakup probability, therefore the surface potential has been 
identified here with the breakup potential. 

Now we discuss the hypothesis leading to Eq.(jD). They have been already 
discussed in Ref. j2] but we report them here too for the sake of completness. 

In this paper we are concerned with reactions performed at energies well 
above the Coulomb barrier where many inelastic channels open at about the 
same distance of closest approach. The effect of the breakup is most im- 
portant at large distances of closest approach {he > Rs), where it represents 
the dominant reaction mechanism. If the breakup probability is needed at 
smaller impact parameters, then the values calculated by perturbation the- 
ory, have to be multiplied by the core survival probability, as discussed in 
Eq.(V.8.1) of Broglia and Winther and also used in relation to halo breakup 
by several authors. The effect of all inelastic channels n different from the 
one we are interested in, can be taken into account by introducing a damp- 
ing factor Pq. Therefore the breakup probability Pb^^ at all distances can be 
defined as 

Pbup = Pbup 11(1 - pn) ^ Pbup exp(- J2 Pn) = PbupPo (5) 

n n 

Each elementary inelastic probability Pn and breakup probability is small 
and pb in particular, can be calculated in time dependent perturbation 
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theory, as done in [221 • reactions with halo projectiles the damping factor 
Pq has also been referred to as the core survival probability after the halo 
breakup or as the core elastic scattering probability. The breakup probability 
Eq.(0) integrated over the core-target impact parameter be has been widely 
used in the literature to get breakup cross sections. 

In this paper we will treat both nuclear and Coulomb breakup as inde- 
pendent process with the formalism used in Ref. namely we will calculate 
nuclear breakup in the eikonal approximation and Coulomb breakup in first 
order perturbation theory. In Ref. 1^ we showed that this is appropriate 
for a one neutron halo since the interference effects are small and the higher 
order effects in Coulomb breakup are negligible. This has been confirmed by 
recent experimental data [HI]. The case of a proton halo or a two-neutron 
halo might need the full all-order approach. 

Then the total breakup probability will be given by 

PKr,ibc)=pUK)+pZ(K). (6) 



2.1 Nuclear breakup 

The optical potential due to nuclear breakup was extensively discussed in 
Ref.[S]. We report here on some of the most important results which will 
help us also constructing the potential for the Coulomb breakup channel. 
The nuclear breakup probability given by the eikonal model is 



<(M ^ ^/^/rfb„(|l-5(b„)r + l-|5(b„)p) 

X (^^El<^Ubn-be,A;Or (7) 

where C'^S is the spectroscopic factor for the initial state, be and b^ are the 
core-target and the neutron-core impact parameters respectively. The calcu- 
lation of the nuclear breakup probability is done here in a reference frame 
with the center at the origin of the target. It is important to remark that 
the above expression takes into account to all orders the neutron target final 
state interaction via an eikonal S-matrix. In this way neutron elastic scat- 
tering and absorption on the target are treated consistently. The integrand 
in Eq.(7) is surface peaked and goes rapidly to zero at the interior of the 
projectile potential (cf. Fig. 5 of Ref. j^). This is because of the natural cuts 



6 



Table 1: Energy dependent optical model parameters for the n-Pb interac- 
tion, ay = 0.55/m, aw = 0.3/m , ry = 1.25/m, = 1.26/m. 



Energy 


V 


Wv 


Ws 


MeV 


MeV 


MeV 


MeV 


20 


-44.8 


-2.03 


-11.02 


40 


-42.8 


-4.02 


-8.50 


80 


-38.8 


-6.08 


-4.87 



introduced by S{bn) and because of the large values of be- Thus the initial 
wave function is needed only at large radii and it can be approximated by 
its asymptotic form which is an Hankel function 

<Pim{r) = -t'Q-fh^\t-fr)YiUO,<P), ir»l, (8) 

where Cj is the asymptotic normalization constant. Its one-dimensional 
Fourier transform reads 



m 



e 



-2r?|bn-bc 



^ C'^— r-Pi{X,), (9) 

2?7|bn - bc| 

where Xi = 1 + 2ki/'y. The divergency in the RHS of Eq.ijH} is compensated 
in Eq.((7j) by the rapid decrease of the terms depending on 5'(6„). 

Eq.((7j) is the neutron breakup probability from a definite single particle 
state of energy Si, momentum 7 = y/—2mei/h, and angular momentum 
/ in the projectile to all possible final continuum state of energy Sf with 
respect to the target and momentum kf = ^'ImSf/h. In our notation ki = 
{sf — Si — ^mv'^) I (hv) and k2 = {sf — Si + |mt>^) / (hv) are the z components 
of the neutron momentum in the initial and final state, respectively, rj'^ = 
kf + 'j'^ = k\ — k^ is the modulus square of the transverse component of 
the neutron momentum. Recoil effects on the neutron are thus taken into 
account. 

For the purpose of this paper the important thing is that in Eqs.(j7J 
the main dependence on the core-target impact parameter is contained 
in the exponential factor Q-'^v\bc-b„\ _ After the 6„ and €f integration, the 
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breakup probability has still an exponential dependence on be for impact 
parameters larger than the strong absorption radius. This has been shown 
already in Fig. 4 of Ref.|Sni and Fig. 1 of Ref.[^. Also Eq. © has a 
maximum in correspondence to the minimum value of rj = j. Therefore the 
6c-dependence of the breakup probabihty p^^^{bc) will be of the exponential 
form p^^^{hc) ~ e^'''^/" with a ^ (27)""'^ where 7 is the decay length of the 
neutron initial state wave function. We now assume at large distances, where 
Po = 1 the same exponential dependence for the absorptive potential due to 
nuclear breakup Wg{r) = e~'^^°'. We assume also, as indicated earher 
on, a straight line parameterization for the trajectory r(t) = be + ft, then 
Eq.(j2)) reads 

/ + OO 
^ W^s{^c.z)dz = --vl^{he). (10) 

The LHS can be approximately evaluated as 



/+OO /■+OO z2 , , I 

W^{K, z)dz = / e-^'^^—c^/'^dz = W(f v^27r6,ae-'=/^ (11) 
'Oo J —00 

where we assumed be » z in the second step. Equating the RHS of Eqs. (fTn|l 
and (jlll) and renaming the distance be as r gives 

^^(r) ^ (12) 

Eq. (fT^ shows explicitly, as already discussed in Ref. (5j, that the long 
range nature of the nuclear breakup potential originates from the large decay 
length of the initial state wave function. For a typical halo separation energy 
of O.SMeV, a = (27) = 3.2/m, while for a 'normal' binding energy of 
lOMeV, a = 0.7 fm as expected. Therefore the parameter a will depend 
mainly on the projectile characteristics and not on the target. 



2.2 Coulomb breakup 

Coulomb breakup can be taken into account as well, following the formalism 
of j33j and the first order perturbation theory probability reads 



Pbur 



AC^SC^ 
M%1 



de, 



rrink 



dVLi. 



zuKi (zu) + izuKo (zu) 



km (k) 

(13) 



8 



The constant Cq = (iiZpZ^e^ /hv is a dimensionless interaction strength and 
w = — ei)bc/hv = ubjv is the adiabaticity parameter. There is part of 
the core recoil which goes to center-of-mass motion. This effect is included in 
the effective charge parameter f3i = nin/mp which is a relevant correction for 
light nuclei. The Coulomb breakup probability is best calculated in the pro- 
jectile reference frame. Thus Ek is here the neutron final energy with respect 
to the core. Obviously after integration over the final energy both nuclear 
and Coulomb probabilities are independent on the reference frame they were 
calculated in. The functions Kq and Ki are modified Bessel functions. 

If the initial state wave function is approximated again by its asymptotic 
form Eq.(jSl), then the general form of the initial state momentum distribution 
is given by the three-dimensional Fourier transform of Eq. (jH]) 

0i„(k) = AttQ ..^^\ .. Yi^Uk) (14) 
7'(/i;^ -|- 7^) 

where k = {kx,ky,kz) is a real vector. 

In the amplitude for Coulomb breakup the Fourier transform of the bound 
state wave function is well approximated by the Fourier transform Eq. ()14p 
of the corresponding Hankel function provided that < 1 and kR < 1 
where R is the radius of the neutron-core potential. These conditions are 
often satisfied in the Coulomb breakup of a halo nucleus. For example for 
^^Be, R ~ 2.5/m, 7 = 0.15/m~^ and the largest final momentum entering 
the calculation is \k\ ~ 0.3/m~^. For this range of final momenta the Fourier 
transform Eq. (fT^ is larger by about 10% than the Fourier transform of the 
Woods-Saxon numerical wave function but it has a very similar shape. The 
assumption of the asymptotic form is expected to overestimate the Coulomb 
excitation by about 20%. The nuclear excitation is better approximated since 
the Fourier transform is taken only in one dimention while |bn — bc| is always 
well outside the neutron-core potential. There are however other incertitudes 
which affect the total probability values. One is the not perfectly known 
spectroscopic factor. The other is the value of the asymptotic normalization 
constant which depends on the geometry used for the neutron-core potential. 
This last incertitude would be present even if the proper Woods-Saxon wave 
functions were used. Also it is important to note that the 2s wave function 
has a node which gives a divergency in its Fourier transform around k = 
0.55/m~^. However for the situations studied in this paper the integrand in 
Eq.()13p has negligible values for such large final momenta. Also using wave 
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functions calculated in a square-well potential it is easy to show that for 
\ei\ < IMeV and up to / = 2 the contribution of the internal part of the wave 
function to the full Fourier transform can be neglected. Here the final state 
has been taken as a plane wave. We have shown in Fig. (8) of Ref . [34j that 
using a proper continuum wave function gives negligible differences when the 
initial state has 1 = 0. 

For a / = initial state which is appropriate for the halo breakup of ^^Be 
for example, and after integration of the angular variables, Eq.()13|) reads 



pfJM « 'ic^S /...|^ {.^Kl (.) + .^Kl (.)) 

(16) 

We show in Appendix A that inserting this result in Eq.© leads to an 
expression consistent with Eq. (12) of Andres et al. ^1]. Those authors have 
developed a semiclassical method to obtain a polarization potential for the 
Coulomb breakup channel at energies close to the barrier. Our formalism 
can then be viewed as a natural extension of the formalism of |14j in the 
high energy, straight line trajectory limit in which the eikonal model of the 
phase shift Eq.© is valid. Also our approach can be considered consistent to 
that of Ref. |2H1 since those authors showed the consistency of their method 
to that of Ref. I 



2.3 Extraction of the imaginary potential 

Eqs.(I2I), © and (jlj) show that the imaginary part of the phase shift and 
the modulus square of the S-matrix which takes into account the breakup 
channel can be simply obtained from the breakup probability and it would 
not be necessary to know the actual form of the optical potential. However to 
understand the physical origin of its characteristics and in order to use elastic 
scattering data for structure studies it is useful to deduce an analytical form of 
it. To this goal we have followed the method outlined in Sec. 2.1. We started 
by calculating the total Coulomb and nuclear breakup probabilities from the 
numerical integration over the neutron final continuum energy of Eqs. (|7j) and 
(fT^. The neutron-target optical potential used here to calculate the neutron 
breakup probabilities is from Ref. (HH] and it has the usual Woods-Saxon form 
for the real and imaginary volume parts and a Woods-Saxon derivative form 
for the surface imaginary part. The parameter values are given in Table 1. 
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Figure 1: Nuclear and Coulomb breakup potentials 



We then fitted the total probabihties to a sum of exponentials. This allows 
a simple extraction of the potential from the relation Eq. ljllj) and also a 
transparent physical interpretation of the large diffusness. 



exp(-6c/«n), 

n 

J^Bn exp (-6c//5„), 



(16) 
(17) 



The corresponding parameters are given in Table 2. The parameters 
and obtained show clearly that the exponential tail of the breakup 
probabilities have a very long range, in particular in the Coulomb breakup 
case. To obtain forms of the potentials valid at all distances, we used then 
Eqs.dni), (fTn|l and (fTTjl . renaming again be with r such that 



~PoJ2An exp{-r/an) J' 
—P^Y^Bn exp(-r//?„) J— — 

2 V V^^PnT 



(18) 

(19) 
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Table 2: Best fit obtained for tlie nuclear and Coulomb halo breakup prob- 
abilities Eqs.(f7|). (fT^ for ^^Be incident on ^°^Pb. The parameters a„ and /3„ 
are in fm, while and i?„ are dimentionless. Incident energies in A.MeV 





Nuclear 


Coulomb 






20 


40 


80 


20 


40 


80 






1860.38 


303.23 


339.45 


9.672 


3.922 


1.734 




A2 


19.73 


7.91 


1.91 


1.558 


0.6013 


0.2549 


B2 










0.1531 


0.0587 


0.0248 


B3 


ai 


1.1924 


1.5334 


1.5356 


2.7273 


2.9438 


3.1046 


Pi 




2.5403 


2.7918 


3.4955 


6.5189 


7.4850 


8.2713 


P2 










15.6863 


19.9322 


24.1546 





Table 3: Optical potential parameters for the bare ^°Be, ^^Be-^°^Pb interac- 

1 /3 1/3 

tion. Radii are calculated with the Ri = ri{Ap + Ar/. )fm convention. 



Energy(A.MeV) 


V 


rv 


ay 


Wv 


rw 


(iw 


20 


-80 


1.020 


0.78 


-66.7 


1.110 


0.39 021 


40 


-70 


0.920 


1.04 


-58.9 


0.890 


0.89 Ha! 


80 


-36 


1.087 


0.78 


-40.0 


1.040 


0.41 Plj 



The core survival probability has been parameterized as 

Po(6J = |^^^|2 = exp(-ln2e[(^=-^^)/'*ol), (20) 

where Qq = O.Gfm and Rg = 1.4(y4p'^ + A^^)fm is the strong absorption 
radius |17j. 

3 Results 

In order to sample the quantitative accuracy of the simple analytical model 
presented above we discuss now some numerical examples. The potentials we 
will discuss derive from the breakup of the 2si/2 state of ^^Be, with separa- 
tion energies 0.5MeV^, asymptotic normalization constant Ci = 0.91/m~^/2 
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Figure 2: Comparison between the results of two calculations according to the eikonal 
model (solid line) and the optical model (dotted line) at 20, 40 and 80 A.MeV for ^"Be 
_l_208p^ with the potential 031, right hand figure, and with the potential 01], left hand 
figure. 

and spectroscopic factor C'^S = 0.77. Cj was obtained from a wave function 
calculated in a Woods-Saxon potential with radius parameter r = 1.25/m, 
diffuseness a = 0.8/m and depth Vq = — 52.4MeV", adjusted to fit the neu- 
tron separation energy. 

As we discussed already in detail in Ref . [S] the diffuseness of the nuclear 
breakup potential is a very important parameter. It depends weakly on 
the incident energy and its large value between 2 and 3 fm is due to the 
large decay constant of the neutron halo wave function and thus to the low 
separation energy. In the case of the Coulomb potential we see that the 
diffuseness /?„ is even larger, and furthermore increasing with energy. This is 
clearly a reflection of the long range, energy dependent effects of the Coulomb 
potential itself and of the core-recoil it causes. 

From the diffuseness parameters values shown in Table 2 we notice that 
at each incident energy the largest diffuseness values are 

f3s^{er'-e,)/hv (21) 

where = 2 3MeV is of the order of the largest continuum energy 

for which there is an appreciable breakup probability. This effect is clearly 
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related to the behavior at very large distances of the Bessel functions Kq i ^ 
in Eq.(|13|). Therefore in the case of the Coulomb breakup potential the 
diffuseness value is related to the adiabaticity parameter and it depends on 
the initial separation energy but also, strongly, on the incident energy. Thus 
this potential is much more dependent from the reaction dynamics than the 
nuclear breakup potential. 



T ' 1 1 1 1 1 1 1 1 1 ' 1 1 1 •- 

20 A.MeV 




2 4 6 8 10 12 14 16 1B 20 

CM (deg) 

Figure 3: Angular distribution for the elastic scattering at 20 A.MeV of ^'^Be and ^^Be 
from 208pb. 

We show now in Fig.l the radial shapes of the potentials calculated for the 
breakup from the 2s state of ^^Be in the interaction with '^^^Pb. Results are 
given for three laboratory incident energies: Ej„c=20 A.MeV, 40 A.MeV and 
80 A.MeV. The solid lines are the Coulomb breakup potentials while the dot- 
ted lines are the nuclear breakup potentials. As discussed above our results 
show that the potentials are energy dependent and the Coulomb breakup 
potential has a much longer range than the nuclear breakup potential. 

The optical potential has one of its most interesting application in the 
calculation of elastic scattering angular distributions. Since according to 
Eqs.(j21) (fnH) our phase shifts for the breakup channels are calculated in 
the eikonal model, it seems appropriate to calculate the angular distributions 
within the same model. Therefore we first show in Fig. 2 a comparison of 



14 



40 A.MeV 




Be+ Pb+Coulomb Breakup \ 
Be+^^^Pb+Nuclear Breakup \ \""" 

Be+^°^Pb+Nuclear+Coulomb Breakup \ '• 



1E-3 I 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 

d CM (deg) 
Figure 4: The same as Fig. 3 but at 40 A.MeV 

calculations made with the optical model code OPTIM jJUI and with our 
modified version of the eikonal code DWEIKO jH] for the reaction at the 
same energies at which potentials have been given in Fig.l. For the lowest 
energy DWEIKO contains a further recoil correction obtained by shifting the 
integration impact parameter, as explained in Ref.jH]. The optical model 
parameters for the volume parts of the bare potential are given in Table 3. 
Since there are no data in the literature for the ^°Be+^°^Pb system at the 
incident energies we are concerned with we have tried to use two modified 
potentials obtained from fits of ^^C+^°^Fh and ^He+^o^Pb Hj data. 
It is clear that for both potentials the eikonal model works very well up 
to 10°. At very high incident energies when an optical model code could 
have difficulties in handling the large number of partial waves necessary for 
heavy ions, an eikonal calculation is much simpler and equally accurate. A 
serius problem for the check of accuracy of theoretical models and related 
numerical applications of the kind of problems discussed in this paper is the 
lack of experimental data and of reliable bare potentials for the core-target 
interaction. Therefore the calculations discussed in the following have to be 
considered as exploratory. 

In Figs. 3, 4 and 5 we show then the eikonal model angular distributions 
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for ^°Be+2°^Pb by the thick sohd hne. The thin sohd hne is for ^^Be+^°^Pb 
with the bare volume potential of Table 3. The dashed and dotted lines in- 
clude the Coulomb breakup and the nuclear breakup potentials respectively, 
while the dot-dashed line includes both. As expected in most of the angu- 
lar range the Coulomb breakup reduces the elastic scattering more than the 
nuclear breakup. 
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Figure 5: The same as Figs. 3 and 4 but at 80 A.MeV 

Finally we discuss another important effect namely how the elastic scat- 
tering total probability changes as a function of the impact parameter or 
angular momentum when there is a large breakup. In Fig.6 we show with 
the solid line the core-target S-matrix Sqt- Because of the incertitudes in 
the bare optical potentials used for the angular distribution calculations we 
prefer here to use for Sct a unique form given by Eq. ipnj) . The nucleus- 
nucleus S-matrix Snn, shown by the dot-dashed line, is calculated then from 
Eq.©, which contains the effect of the halo breakup. At a fixed impact pa- 
rameter the effect of the breakup is to reduce the elastic probability given by 
the modulus square of the S-matrix. The unitary limit is attained at much 
larger 6c-values and the reaction cross section receives significant contribu- 
tions from a large range of impact parameters. The reduction is obviously 
more pronounced at the impact parameters larger than the strong absorption 
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Figure 6: Scattering matrix with and without breakup potentials contributions 



radius. The value of the strong absorption radius, defined as \Snn{Rs)\'^ = | 
changes appreciably, increasing of about 1 fm at 20 A.MeV, 0.8 fm at 40 
A.MeV and 0.6 fm at 80 A.MeV. The large increase is mainly due to the 
Coulomb breakup effect. This can be seen by looking at the dotted and 
dashed lines which represent Snn when only the nuclear or the Coulomb 
breakup effect is included, respectively. 

Another significant effect of the imaginary surface potential is seen in the 
reaction cross sections obtained from the S^n and given in Table 4. We 
show also a number of other cross section values necessary to justify the 
phenomenological inputs of our calculations and their consistency. The first 
two rows give the total cross sections for the systems ^°Be + ^osp^ and ^^Be 
_l_ 208p]^ respectively, calculated with the bare optical potentials of Table 
3. The cross sections in the next three rows contain also the effect of the 
nuclear breakup, of the Coulomb breakup and of both at the same time. All 
these cross sections have been calculated in the eikonal approximation with 
the code DWEIKO. There is an increase of the order of 500-7-600 mb with 
respect to the bare (no breakup) optical potential when the nuclear breakup 
potential is taken into account. When the Coulomb breakup is included we 
find an increase of about 1.6 b to 4.4 b depending on the incident energy. 
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The total increase in the reaction cross section is very large, corresponding 
to a variation of a factor from 2 to 3 at some energies. Similar increases in 
the cross sections have been found in Ref.p^^. These values are consistent 
with the nuclear and Coulomb breakup cross sections obtained from the 
integration over the core-target impact parameter of the probabilities Eqs. (|7j) 
and (jl3j) including the core survival probability Eg . (1201) and given in the 
next two rows. Finally the free particle n-Pb elastic, inelastic and total cross 
sections obtained in the eikonal model by using the potentials of Table 1 are 
given. The last row contains the experimental free particle cross sections. 
As expected at 20MeV our free particle cross sections are still not very close 
to the experimental values, but already at 40MeV the eikonal approximation 
works very well. 

As discussed already in the relation between the total reaction cross 
section and the breakup cross section, can be understood by expanding the 
exponential in Eq.(0} to first order in P^^^ and integrating over the impact 
parameter he- One immediately finds 

= l-\ScT{hc)? + 

+ \ScT{hc)\\plS^c)+pl,{K)) (22) 

ctnn = 2n lbJK{l-\SNNibc)\^) (23) 

^ acr + < + (24) 

The fact that the cross sections values of Table 4 are in very good agree- 
ment with the relations contained in Eqs. (j23|) and (j24j) is a proof of the 
accuracy of the hypothesis Eq.(jD) for the nucleus-nucleus S-matrix and of 
the separation of the imaginary potential into a volume and a surface term, 
with the surface term identified with the halo breakup potential. 

4 Conclusions 

In conclusion we have presented a simple analytical method to obtain the 
surface component of the imaginary part of the nucleus-nucleus optical po- 
tential to be used in the elastic scattering calculations between a halo or 
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Table 4: Cross section values (in mb) discussed in this paper. See text for 
details. 





Energy(A.MeV) 


20 


40 


80 


^^Be +^°* Pb total cross section 


CTCN 


2362 


2897 


2537 


Be +208 p5 total cross section with bare potential 


O'NN 


2644 


2969 


2593 


(Jnjs[ including nuclear breakup 




3185 


3471 


3060 


(Tnn including coulomb breakup 




6685 


5477 


4152 


aNN including coulomb and nuclear breakup 




7074 


5900 


4584 


nuclear halo breakup 




615 


551 


490 


Coulomb halo breakup 




4280 


2638 


1570 


jj_208pj, elastic free particle 




2554 


2623 


2467 


n-2"®P6 inelastic free particle 




2018 


1902 


2097 


jj_208p^ total free particle 


^Tot 


4548 


4437 


4867 


j^_208p^ experimental 


O-exp [455 


5888 


4392 


4817 



weakly bound nucleus and a heavy target. The surface potential is due to 
Coulomb and nuclear breakup. The main purpose here was to relate the char- 
acteristics of the potential to the special properties of the breakup channels 
for weakly bound nuclei. At high incident energy (> 20 A. Me V) the evalua- 
tion of the potentials amounts in fact just to the calculation of the breakup 
probability as already shown in Ref . |Sj . If breakup from core excited states is 
to be included, then it suffices to sum up the relative probabilities according 
to Eq.©. 

The method to include Coulomb breakup is an extension of that pre- 
viously used to calculate microscopically the effect of transfer and nuclear 
breakup channels on the imaginary potential. The shape of the surface imag- 
inary potential and its parameters are determined univocally by the shape 
of the breakup probability. An interesting result is that the diffuseness of 
the potential for the nuclear breakup reflects the decay length of the valence 
neutron wave function and therefore it depends mainly upon the projectile 
structure, but not so much on the reaction dynamics. For the Coulomb 
breakup there is instead a strong dependence on the dynamics which leads 
to an increase of the diffuseness value when increasing the incident energy. 
The strength is also energy dependent. Analytical calculations contained in 
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appendix A have shown that the potential proposed here is consistent with 
and can be viewed as an extension to high energy, of other theoretical models 
developed at energies close to the Coulomb barrier The numerical cal- 
culations on the other hand show consistency with the works of other authors 
[T3] for similar, light halo systems. Furthermore we have given an explicit 
justification for the long range of the polarization potential. 

From Feshbach 50. formalism it is known that related to the imaginary 
potential, which comes from the second order, complex term, there is also 
a real correction to be added to the first order term, the folding potential. 
However it has been shown by the authors of Ref.^J] that the Coulomb 
excitation polarization potential is purely imaginary at high energy. On 
the other hand for the nuclear breakup potential, we have shown in Ref.jH] 
that the second order real correction is small and negligible, while other 
calculations [HIE] for ^^Li projectile found a not so small real polarization 
potential. To make sure that the second order real correction is negligible 
also for heavy targets such as ^^'^Pb, we have added phenomenologically a real 
repulsive part to our microscopically calculated imaginary potential, with the 
same exponential form, same diffuseness and variable strength, but we found 
no noticeable differences in the angular distributions. Two open questions 
remain then to be addressed for future work. One is which observable will 
be more sensitive to the real correction term. The other is a systematic 
microscopic calculation of such a term for a number of different systems and 
energies. 
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A Connection to low energy approaches 

In our formalism the imaginary part of the optical potential due to Coulomb 
breakup is given by 

/+CO fi 
^ Wsm)dt = --pl^, (25) 

while Eq. (12) of Andres et al. [Hj reads 

/+00 
Ws{r{t))dt = ha. (26) 
-oo 

Where both and a depend on the incident energy and on the classical 
trajectory of relative motion. In our case such a trajectory is a straight 
line since our approach applies to incident energies well above the Coulomb 
barrier. 

The two approaches are consistent if a ^ — p^^/2 in the high energy 
limit. To show that this is true we write explicitly the expression in Ref. ^1] 

2 



a 



where are the well known Coulomb integrals jlo] which can be expressed 
in terms of Bessel functions of imaginary order. However according to Eq.(28) 
of Ref. ^48j in the high energy limit 

I{E1, ±1) = Ji±i = ^zuK,{zu), (28) 



where ac is the Coulomb length parameter and now Ki is an ordinary Bessel 
function of real index. 

On the other hand considering Eq. ljlSp for we remark that it is well 
known and shown for example in Fig.l of Ref.|ini that zu'^Kq {w) is much 
smaller than vu'^Kf (cij) for values of w ^ 0.1, which happens for heavy ions 
at high energies. We can then write Eq.()15p as 

since for the case studied in this paper the explicit form of B(E1) is 

de, - ^ ^^(C'^/^i^^^) -2j^TW ^ ^ 
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which is the consistent with Eq.(6.5) of Ref.|m] or Eq. (7.24) of [2H1 and 
it is the exphcit form of the B(E1) obtained for an s-initial state using the 
asymptotic form of the wave function. 

Finally using Eq. ()28|) in Eq. ljTTj) and comparing with Eq. ()29|) we obtain 
that a = — p^^/2 in the high energy, straight line trajectory limit. However 
we remark that our probability Eq. (jl3|) is more accurate than the high energy 
limit of a of Eq. (j27j) since it contains also the term representing longitudinal 
excitations, proportional to Kq and to kz, the parallel component of neutron 
momentum. 
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